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ABSTRACT 


A  model  of  the  water  gain  of  empty  structure!  stored  In  the  open  is 
derived.  Using  this  model  and  a  modified  Monte  Carlo  procedure,  a 
statistical  analysis  is  made  of  probable  water  gains  at  12  stations  by 
enclosures  which  are  free  breathing  and  which  have  check  valves  set  at 
+2  -1,  +1  -1,  +1  -1/2  and  +1/2  -1/2  psig. 

It  is  concluded  that  a  breather  charge  of  one  ounce  per  cubic  foot 
should  give  about  six  months  storage  life  anywhere  in  the  world,  as¬ 
suming  adequate  container  drawdown.  Even  the  +1/2  -1/2  psig  valve 
required  no  additional  desiccant  although  it  might  be  expected  to  breathe 
as  much  as  38  times  per  year  at  Las  Vegas,  Nevada  with  2d  confidence. 

The  least  pressure-vacuum  spread  which  will  not  adversely  affect 
a  reasonable  storage  life  was  not  discovered. 
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INTRODUCTION 


A  shipping  container  in  the  open  is  exposed  to  deteriorating 
inAuences  such  as  sunlight,  precipitation,  water  vapor,  sand  and 
dust,  microorganisms,  etc.  Design  criteria  for  the  various  extreme 
conditions  have  been  published.  *  U) 

So  long  as  protection  of  the  contents  depends  upon  &  sealed 
barrier,  extreme  condition  values  are  useful  design  criteria.  When, 
however,  the  container  is  deliberately  designed  to  inhale  and  exhale 
air,  reliance  on  extreme  conditions  produces  unrealistic  designs. 

The  gas  laws  teach  that,  in  a  completely  sealed  container, 
the  pressure  will  vary  with  the  internal  temperature.  In  designing 
such  &  container  to  be  able  to  go  anywhere  in  the  world,  it  is  obvi¬ 
ously  sufficient  to  design  for  the  estimated  maximum  temperature 
which  would  be  encountered.  The  estimated  minimum  low  tempera¬ 
ture  will  dictate  the  maximum  vacuum. 

For  a  breathing  container,  however,  maximum  pressure  is 
secondary.  **  The  mass  of  air  exchanged  with  the  surroundings  is 
the  governing  factor  since  dry  air  exhaled  will,  as  temperature  drops, 


*  See  literature  references  at  end  of  this  report. 

**  On  controlled  breathing  containers,  the  relief  valve  settings  will 
determine  the  pressures  and  vacuums  developed.  On  free 
breathers,  back  pressure  in  tubing  or  desiccant  beds  will  be  con¬ 
trolling. 
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be  replaced  with  air  whose  water  content  is  characteristic  of  the 
environment.  Hence  it  is  necessary  to  evaluate  what  the  internal 
temperatures  might  reasonably  be  expected  to  be  in  typical  cli¬ 
mates.  With  such  internal  temperature  data  it  is  possible  to: 

a.  Determine  desirable  valve  settings  for  controlled 
breathing  containers. 

b.  Evaluate  the  total  water  intake  load  which  must  be 
taken  up  by  desiccant  in  the  container  and/or  in  the  breathe r* 

This  study  was  undertaken  to  obtain  these  design  criteria. 
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METHOD  OF  COMPUTING  INTERNAL  TEMPERATURES 


In  order  to  standardize  results,  this  report  assumes  that  one 
vertical  surface  of  the  container  is  exposed  normal  to  the  solar  azi¬ 
muth  at  time  of  equilibrium. 

For  thin  walled  structures,  such  as  containers,  Parmelee  and 
Aubele  (2)  8how  that  instantaneous  heat  transfer  rate  (q/a  =  heat 
per  unit  area  per  unit  time)  on  the  sun  side  can  be  computed  from: 

•J-  -  Uj  *“t|  )  (1) 

where 

Uj  =  ox'er-all  transmission  coefficient  on  the  sun  side 
te  =  the  sol-air  temperature  on  the  sun  side 
t}  =  the  inside  air  temperature 

(3) 

The  sol-air  temperature  is  v  '  that  temperature  which  would 
exist  if  the  entire  heat  transfer  process  werts  occurring  by  conduction. 

Similarly,  on  the  shaded  side,  which  is  losing  heat,  the  heat 
transfer  function  is 

I  *  u,  Oi  -t„h>  <2> 

where 

U2  =  over-all  transmission  coefficient  on  the  shade  side 
*esh  =  the  sol-air  temperature  in  the  shade 
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In  solving  specific  problems,  it  is  necessary  to  consider  known 
areas.  For  purposes  of  this  study  we  consider  a  rectangular  cross- 
section  container  whose  dimensions,  reasonably  typical  of  many  designs, 
are  L  =  2W  =  2H. 

With  such  a  container  the  area  of  the  top  is  the  same  as  the  area 
of  the  side  and  the  combined  area  of  the  ends  equals  the  area  of  the  top. 
For  simplicity,  exchange  between  the  container  and  the  ground  is  ig¬ 
nored.  This  is  usually  an  outward  flow  during  the  day.  Thus,  this 
omission  will  result  in  computing  slightly  higher  internal  temperatures, 
thus  producing  a  conservative  effect  on  water  load. 

Under  these  assumptions,  the  total  area  losing  heat  equals  the 
total  area  gaining  heat.  Making  the  further  assumptions  that  equilib¬ 
rium  is  reached,  *  and  that  one  vertical  wall  of  the  container  is  normal 
to  the  sun's  azimuth,  we  may  write: 


ts  = 


t 


u, 


e*h 


I  +■ 


u, 

ua 


(3) 


*  This  assumption  is  reasonable  for  empty  containers.  It  ignores, 
however,  the  heat  capacity  of  container  contents.  For  precise  com¬ 
putation  in  specific  instances,  however,  one  could  make  the  com¬ 
plete  analysis  using' Fourier  integrals.  The  effect,  here,  is  judged 
conservative. 
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4.  /*  f  *  "  S'  ^  ^ 


Now,  Parmelee  and  Aubele  ^  define  t# 


t.  =  *ito-llS!_  +  t0  =  At#+  t( 


where 


>  solar  absorptivity 

=  emissivity  at  surface  temperature  T§ 

=  average  direct  solar  radiation  rate  on  the  top  and  the 
one  exposed  side. 

=  radiation  exchange  factor  on  sun  side,  made  up  of  the 
balance  between  incoming  sky  radiation  and  reradiation 


=  convection  coefficient  on  sun  side 


a  radiation  coefficient  on  sun  side 


Similarly, 


t«.h  -  ♦ 


«R« 

^S0f+  hR» 


-  t  ~  Ati#h 


where  the  subscript  2's  refer  to  shade  side  coefficients. 
Combining  (3),  (4)  and  {5),  we  have 


=  t«  + 


'*£ 


Values  of  Uj,  hCQ  and  hy 

The  over-all  transmission  coefficient  is  a  somewhat  variable 
number  made  up  from  convection  and  radiation  coefficients  on  both 
sides  of  the  wall  and  reflecting  the  thermal  capacity  of  the  wall  ma¬ 
terial.  For  the  sun  side  we  have 


and  for  the  shade  side 


(7) 


where 


(8) 


=  inside  surface  coefficient,  actually  composed  of  values 
representing  still  air  convection  and  radiant  exchange 
with  the  contents. 


d  =  wall  thickness 


k  =  thermal  conductivity  of  wall 

h  =  outside  convection  coefficient 

co 

h  =  outside  radiation  coefficient 

r 


Inasmuch  as  the  contents  are  unknown  in  this  case,  it  is  suf¬ 
ficient  for  engineering  purposes  to  use  the  generally  accepted  value  h| 
1.65  which  reflects  average  conditions. 


8 


In  thin  walled  aluminum  alloy,  steel  or  even  plastic  containers, 
it  is  readily  verifiable  that  the  ratio  d/k  is  usually  considerably  less 
than  .  003  and,  therefore,  can  be  ignored.  * 

The  value  selected  for  hco  will  vary  with  wind  velocity,  all 
other  things  being  equal.  For  purposes  of  this  study  still  air  is  as¬ 
sumed  to  prevail.  This  condition  was  selected  as  being  the  most  severe. 
The  average  film  coefficient  for  the  sun  surfaces  was  taken  ^  as 


'coi 


0.33 


<T„  -  V 


o.ts 


(9) 


On  the  shade  si^e,  only  vertical  surfaces  are  involved  so  that 


(10) 


where  ts  represents  the  surface  temperature. 

Few  data  are  available  on  surface  temperatures  of  objects  ex¬ 
posed  to  direct  sunlight.  There  are  some  indications  that  they  can  be 
very  high.  For  example,  ts  =  198°F.  when  tQ  =  120°  have  been  re¬ 
ported  from  measurements  of  aircraft  wing  surfaces  at  Tucson.  ^ 

In  order  to  evaluate  hcQ,  tg^  was  assumed  to  be  65  percent  higher  than 
tQ.  Shadey  side  data  are  lacking  but,  on  thin  walled  structures,  they 
must  be  somewhat  higher  than  ambient.  We  have  assumed  t a  constant 


* 


For  example,  low  carbon  steel  has  a  k  of  approximately  26.2  BTU/ 
(hr. )  (ft.  2)  (deg.  F.  /ft. ).  For  1/4"  thickness,  d/kj  =  0.0008. 
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fiv9  percent  higher  than  t0. 

In  computing  sol-air  temperatures,  Parmelee  and  Aubele^ 
define  the  radiation  coefficients  by 


=  «<r(T„4-T04)  t  .15916  (T,,4—  T0 4  )  _ 


T  -  T 

'si  'o 


Tt,  -  T„ 


«a(T„4-T04)  ,15916  (V-T,4) 


T»-  T0 


T..  “  To 


where 


<r  =  Stefan-Boltzman  constant,  0.  173  x  10"®  BTU/hr. /ft^/  °R 
T  =  temperature,  °R 


For  many  paints  used  on  specialty  containers,  and  for  corrugated 
fibreboard  containers,  low  temperature  emissivity  is  approximately  the 
same,  i.  e.  e  =  0.  92. 

Values  of  Uj  and  U£  were  then  determined  graphically  and  their 
ratios  were  also  determined  graphically.  These  ratios  were  found  to  be 
essentially  constant  over  a  wide  range  of  temperatures  with  the  following 
numerical  values: 


=  0.83 
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Ul 

U2 


1.2) 


Further  Simplification  of  tj  Estimating  Equation 

The  symbol  A  te  defined  in  equation  (4),  contains  QC  ( 

the  solar  absorptivity.  For  most  surfaces,  the  absorptivity  and  emis- 

* 

sivity  are,  to  all  intents  and  purposes,  the  same.  Letting  CC  m  c  *0.92, 
permits  us  to  write 


+  (Idov  ”  Ri)“  R4 


(13) 


where 


* 


0.6027 


4> 

'i 


_  0.4163 
^coa  +  ^Ra 


Values  of  <t> 
to  compute  the  various 


and  <j)  are  plotted  in  figure  1  and  were  used 
average  t^. 


*  The  significant  exceptions  are  white  paint,  aluminum  pigmented 
paint  and  unpainted  corrosion  free  surfaces. 
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Computing  R  j  and 


Brunt  (6)  estimates  incoming  scatter  radiation  R8  by 

R8  *  «tT04(q  +  byp#"Ht )  (H) 

where 

a,  b  =  empirically  determined  constants  tabulated  in  Table  I. 
P8  =  saturation  vapor  pressure  at  T0. 

Hj,  s  decimal  relative  humidity  at  T0. 

TABLE  I 

VALUES  OF  CONSTANTS  a  and  b 
(After  Dines7) 


a 

b 

Vertical  surfaces 

0.30 

0.  165 

Horizontal  Surfaces 

0..50 

0.34 

Average  for  Equal  Area  Surfaces 

0.40 

0.253 

The  net  reradiation  factor,  R^,  is  the  difference  between  black 
body  radiation  at  TQ  and  Rg,  or 

R.=  aTo4(l-a-bv^“)  (15) 

Hence,  we  have,  on  the  sun  side 

R(  =  T04(.I04-  .043769v^h*2  )  *  10  (16) 
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and  on  the  shade  side 


R2  *  T04(.I2I  -.02  845 -/pT^)  *IO-# 


(17) 


Equation  (15),  and  its  explicit  forms  (16)  and  (17),  are  tedious 
to  compute  by  hand  for  a  large  number  of  TQ  and  H.  Note,  however, 
that  pg  is  also  a  function  of  TQ  and  we  may  write,  in  functional  form 


R|  =  M  “/Wj 


(18) 


in  which 


f,  (t0)  *  .104  (T0)4  X  I0*# 

f2  (t0)  =  .121  (T0)4  X  io- • 
f8  (t0)  *  .0438  (T0)47p“  XI0-# 

f4  (to)  5  0285  (T0)4ypT  *  I0-* 


Values  for  f^,  f^  and  f^  are  plotted  in  Figure  2. 

Ho  and  t  are  taken  from  the  climatic  data,  p  from  the  steam 

“  O  8 

tables,  and,  with  the  aid  of  figure  2,  Rj  and  R2  are  estimated  quickly. 


Computing  Incident  Direct  Radiation 

With  the  vertical  wall  normal  to  the  sun's  azimuth 

^  *  !dn  C0S  3  1 

\i  *  X*  Si"  *  * 


(19) 
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and,  for  equal  area  horizontal  and  vertical  surfaces, 


I 


dav 


sin#  +  cos#  . 

2  *dn 


(20) 


whore 

=  direct  radiation  rate  on  horizontal  surface 

*2d  =  direct  radiation  rate  on  vertical  surface 

1^  =  direct  radiation  rate  on  a  surface  normal  to  the  sun's 

rays 

ft  =  sun* 8  altitude  at  the  time  and  place  under  consideration. 


Actual  direct  solar  radiation  (I  £jn)  reaching  a  surface  normal  to 
the  sun's  rays  is  a  highly  variable  quantity  whose  value  is  influenced  by 
the  length  of  the  atmosphere  traversed  by  the  sun's  rays  and  the  com¬ 
ponents  of  that  atmosphere,  particularly  water  vapor,  dust,  and  carbon 
dioxide,  industrial  haze  and  the  amount  of  cloud  cover. 

In  this  study  a  clear  day,  free  of  industrial  haze  is  assumed. 
Thus,  the  reducing  variables  of  cloud  cover  and  industrial  haze  become 
conservatively  eliminated. 

Water  content  of  the  atmosphere  at  any  given  time  is  not  reftdily 
predictable.  Moisture  distribution  in  the  atmosphere  is  only  approxi¬ 
mately  measured  by  ground  conditions.  Moon  (8)  calculated  radiation 
reaching  the  surface  using  20-mm  of  precipi table  water  and  estimates 
that  error  would  be  small  for  other  values  of  water  content  in  our  calcu¬ 
lations  . 
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Using  the  foregoing  assumptions,  direct  normal  incident  radiation 
can  be  computed  by 

Idn=TIos429T  (ZD 

where 

T  =  atmospheric  transmission  coefficient 

IQ  =  normal  radiation  at  the  top  of  the  atmosphere.  While 
somewhat  variable  with  the  radius  vector,  the  usual 
value  assigned  is  429  BTU/ft.^/hr. 

Values  of  T  have  been  computed  by  Kimball^)  and  are  repro¬ 
duced  in  the  Smithsonian  Tables  (l0),  For  hand  work,  the  latter  ref¬ 
erence  is  preferred.  For  computer  work,  however,  it  was  found  that 
the  non-linear  T  function  could  be  matched,  in  the  range  of  optical  air 
masses  to  be  considered, by  the  linear  function 

T  =  0.86  -  0.09  m  (22) 


v,  lie  re 

m  =  optical  air  mass 
with  error  less  than  0.01. 

For  sun  altitude  greater  than  10°  (  (3  >  10°),  the  optical  air 

mass  is  computed  from 

m  a  p/D  csc/3  (23) 

Ko 

* 
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where 

p/po  =  station  altitude  correction  factor  taken  from  the  NACA 
atmosphere  tables  02) , 

Combining  (20)  through  (23),  we  have 

Idfly  =  214.5  (sin£  +  cos£)(o.86--i£2 - j  (24) 

p/„  sin  /3 

Note,  here,  that  negative  answers  are  trivial  and  should  be  re¬ 
corded  as  zero.  So  long  as  ft  i  6°,  answers  will  be  positive. 

Solar  altitude,  ft  ,  can  be  read  directly  from  tables  (13),  Qr 
in  a  computer,  by  solving  the  well  known  relationship 

sin  £  csin  L  sin  8  +  cos  L  cos  8  cos  h  (25) 

in  which 

L  =  station  latitude 

8  =  sun's  declination 

h  =  hour  angle  at  time  of  equilibrium. 

Now,  limited  studies  (14)  indicate  that  maximizes  around  3:00  P.  M. 

local  sun  time.  Hence,  for  this  comparative  study,  it  is  sufficient  to 

use  cos  h  =  0.  7071.  In  using  equation  (25),  particularly  in  a  computer, 

it  is  well  to  remember  t(hat  8  can  have  negative  values  to  the  observer 

when  the  sun  is  on  the  opposite  side  of  the  equator. 

The  value  of  8  can  be  estimated  from  tables  of  the  ephemeris 
(15) 

of  the  sun  .  As  might  be  expected,  it  is  a  sine  function  with  half- 
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period  extending  from  equinox-to-equinox  and  with  peak  amplitude 
jf23.5°  (0.  4093  radians),  occurring  at  the  solstices.  Thus,  the  abso¬ 
lute  value  of  S  on  the  i  th  day  (d[)  after  the  spring  equinox  can 
be  computed  from 


was  computed  for  Lincoln,  Nebraska  with  input  values  of  average  high 
tQ  =  100°F.  and  H  =  0.48  for  the  month  of  July  taken  from  the 
world  wide  climatic  summaries  U6).  sky  radiation  on  a  horizontal 
surface  was  found  to  be  143  BTU/sq.  ft.  /hr.  The  value  for  for 
July  17  using  the  summarized  procedure  is  192  BTU/sq.  ft.  /hr.  Thus 
the  total  horizontal  radiation  is  335  BTU/sq.  ft. /hr.  Mackey  for 
the  same  location  and  month,  reports  measured  values  of  337  BTU/ft.  / 
hr.  at  2:00  P.M.,  C.S.  T.  and  293  BTU/ft.  2/hr.  at  3:00  P.M.,  C.  S.  T. 
Thus  it  is  apparent  that  the  computed  values  are  as  close  as  can  be 
expected  to  observed  values  and  that  error  is  on  the  conservative  side. 

It  should  be  emphasized  that  I^av  is  not  necessarily  the  maxi¬ 
mum  direct  radiation  rate  but  is  rather,  the  radiation  intensity  at  about 
the  time  of  maximum  t^  and  tQ.  tQ  always  tends  to  lag  peak  insolation 
(occurring  at  noon  local  solar  time)  because  most  of  the  heat  of  the  at¬ 
mosphere  is  the  result  of  reradiation  in  the  longer  wave  lengths  from 
the  earth's  surface. 
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TABLE  II 


DISTRIBUTION  OF  STATIONS  ANALYZED 


Climate 

Class 

Description 

No. 

of 

Stations 

Climate 

Class 

Description 

No. 

of 

Stations 

Am 

Monsoon 

7 

C8 

Dry  Summer  Subtropical 

19 

A£ 

Tropical  Rain  Forest 

13 

Da 

Humid  Continental,  warm 

summer 

5 

Aw 

Tropical  Savanna 

14 

H 

Undifferentiated  highland 

2 

Bw 

Desert 

12 

Cb 

Marine,  cool  summer 

9 

Bs 

Steppe 

12 

Db 

Humid  Continental,  cool 

summer 

3 

Ca 

Humid  Subtropical 

19 

Dc 

Sub  Arctic 

5 

w  *  j 


SURVEY  OF  INTERNAL  TEMPERATURES 


Selecting  Locations  to  be  Analyzed 

To  solve  the  problem  at  hand,  it  is  necessary  that  locations 
analyzed: 

A.  Be  reasonably  representative  of  the  various  world  climates. 

B.  Be  a  reasonably  possible  storage  site  (e.  g.  the  top  of  Mt. 
Washington  is  not  considered). 

C.  Have  adequate  weather  data  with  statistics  extending  over 
several  years. 

Trewartha's  (18)  modification  of  the  Koppen  climate  classifi¬ 
cation  system  is  the  classif  cation  used  in  this  paper.  Frequency  dis¬ 
tribution  of  the  120  stations  analyzed  with  respect  to  12  climates  is 
shown  in  table  II.  Within  the  climate  types,  stations  were  selected  with 
a  view  to  obtaining  world  wide  distribution,  insofar  as  practicable. 

Basic  source  of  weather  data  for  the  stations  selected  was  the 
world*  wide  tables  (16),  These  data  are  essentially  official  weather 
station  data  and  do  not  necessarily  represent  the  micro-climates  close 
to  the  ground  in  proba'ble  storage  locations.  They  are,  however,  the 
only  readily  available  world  wide  compilation  and,  in  the  absence  of 
more  detailed  data,  must  be  accepted  as  representative  of  true  condi¬ 
tions. 

As  a  preliminary  step,  monthly  average  internal  temperatures 
were  determined  for  these  stations  using  the  middle  of  the  month  for  $  . 


These  results  were  used  to  compute  "average"  water  gain  at  each 
station  for  free  breathing  containers  in  order  to  characterize  the  dis¬ 
tribution  of  these  water  gains  by  climate  types. 

Detailed  computing  results  for  with  intermediate  answers, 
for  one  station  are  shown  in  table  III.  Table  IV  contains  average  t^ 
for  all  stations  analyzed. 
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TABLE  III 


DETAILED  TEMPERATURE  COMPUTATIONS  FOR  OKLAHOMA  CITY 
Latitude  35°  Altitude  1254' 


Property 

May 

June 

July 

August 

Sept. 

Oct. 

Average  Ambient  Hi  (t  ) 

o 

78 

87 

92 

92 

85 

73 

Average  R.H.  at  to 

.57 

.55 

.48 

.46 

.51 

.55 

Sun's  Declination 

19°  9' 

23°  22' 

21°  21* 

13°  41* 

2°34' 

-8°  58' 

Sun’s  Altitude 

47 

49 

49 

44 

37 

29 

Atmospheric 

Transmissibility 

.72 

.72 

.72 

.71 

.69 

.65 

Optical  Air  Mass 

1.30 

1.26 

1.26 

1.37 

1.58 

1.96 

Normal  Radiation 

309 

309 

309 

305 

296 

279 

Average  Direct 

Radiation  (I .  ) 

dav 

218 

218 

218 

216 

207 

190 

0i 

.249 

.237 

.231 

.231 

.240 

.252 

0a 

.301 

.288 

.281 

.281 

.291 

.301 

Ri 

54 

53 

52 

52 

53 

56 

Ra 

80 

82 

83 

83 

82 

79 

Temperature  Gain 

41 

39 

38 

38 

37 

34 

Temperature  Loss 

24 

24 

23 

23 

24 

24 

Average  Max.  Int,  Temp  t. 

95 

102 

107 

107 

98 

83 

TABLE  IV 

AVERAGE  INTERNAL  HIGH  TEMPERATURES  -  6  HOTTEST  MONTHS 


St  a 

No 

Name 

BiH 

2 

3 

t. 

1  4 

5 

6 

1 

Belem,  Brazil 

16 

Am  Climates 
July  103 

104 

105 

105 

105 

104 

2 

Cochin,  Indochina 

43 

Dec 

100 

100 

102 

104 

105 

104 

3 

Colombo,  Ceylon 

25 

Jan 

98 

100 

102 

102 

101 

100 

4 

Georgetown,  B.Gui, 

29 

June 

101 

101 

102 

103 

103 

102 

5 

Manila,  P.I. 

16 

Mar 

105 

108 

108 

106 

102 

103 

6 

Monrovia,  Liberia 

3 

Dec 

100 

101 

101 

103 

103 

103 

7 

Rangoon,  Burma 

60 

Dec 

96 

98 

103 

109 

106 

109 

8 

Belize,  Br.  Hond. 

27 

Af  Climates 
May  103 

104 

104 

104 

102 

101 

9 

Kieta,  Solomon  Is. 

9 

Oct 

102 

101 

102 

101 

102 

102 

10 

Kuching,  Borneo 

5 

Apr 

90 

90 

91 

90 

91 

89 

11 

Menado,  Celebes  Is. 

12 

May 

102 

102 

101 

103 

105 

105 

12 

Noumea,  New  Caled. 

22 

Nov 

97 

101 

101 

100 

101 

95 

13 

Penang,  Malaya 

48 

Feb 

104 

105 

105 

104 

103 

103 

14 

Rio  de  Janiero,  Br. 

26 

Nov 

96 

98 

100 

100 

97 

93 

15 

San  Juan,  P.R. 

43 

June 

100 

100 

100 

101 

100 

99 

16 

Stanleyville,  Congo 

3 

Jan 

104 

104 

105 

104 

103 

102 

17 

Surigao,  P.I. 

15 

Apr 

103 

104 

104 

104 

104 

104 

18 

Suva,  Fiji  Is. 

31 

Dec 

101 

102 

101 

100 

96 

92 

19 

Tamatave,  Madag, 

20 

Nov 

95 

96 

97 

97 

96 

96 

20 

Tarawa,  Gilbert  Is. 

5 

July 

101 

103 

103 

104 

103 

101 

21 

Balboa  Hts,,  Panama  33 

Aw  Climates 
Dec  100 

101 

103 

105 

105 

103 

22 

Bangkok,  Siam 

16 

Feb 

105 

108 

110 

108 

106 

106 

23 

Bombay,  India 

60 

May 

105 

104 

101 

100 

99 

103 

24 

Brazzaville,  F.Eq.Af 

.  7 

Dec 

103 

104 

105 

106 

106 

104 

25 

Calcutta,  India 

60 

May 

99 

104 

102 

101 

99 

98 

25  - 


TABLE  IV  (CONT.) 

AVERAGE  INTERNAL  HIGH  TEMPERATURES  -  6  HOTTEST  MONTHS 


5ti“ 

No__ 

Name 

Data 

"""  let*1 
Mo 

.  1 

2 

t 

3 

!‘  4 

5 

6 

26 

Honolulu,  Hawaii 

26 

June 

96 

98 

98 

97 

93 

89 

27 

Madras,  India 

60 

May 

108 

112 

111 

108 

108 

107 

28 

Manaus,  Brazil 

8 

Jul 

103 

106 

107 

107 

106 

105 

29 

Mandalay,  Burma 

20 

Mar 

99 

106 

105 

101 

101 

101 

30 

Mombasa,  Kenya 

45 

Nov 

96 

96 

97 

98 

98 

96 

31 

Nairobi,  Kenya 

12 

Nov 

93 

93 

94 

96 

96 

94 

32 

Natal,  Brazil 

6 

Dec 

100 

102 

101 

101 

100 

100 

33 

Port  Darwin,  Aus, 

53 

Oct 

110 

107 

105 

104 

'  04 

106 

34 

Saigon,  Indochina 

11 

Jan 

100 

101 

104 

107 

105 

104 

Bw  Climates 

35 

Bahrein,  Pers.Gulf 

16 

May 

106 

109 

112 

113 

108 

100 

36 

Eilat,  Israel 

4 

May 

105 

109 

110 

112 

107 

98 

37 

Isfahan,  Persia 

22 

May 

98 

103 

108 

106 

100 

87 

38 

Karachi,  Pakistan 

43 

May 

109 

108 

107 

104 

103 

102 

39 

Kazalinsk,USSR 

11 

May 

?U 

99 

102 

97 

84 

63 

40 

Kuwait  City  .Arabia 

12 

May 

107 

110 

111 

114 

108 

98 

41 

Las  Vegas,  Nevada 

10 

Apr 

92 

101 

106 

111 

111 

102 

42 

Lima,  Peru 

15 

Nov 

92 

95 

98 

100 

100 

96 

43 

Phoenix,  Arizona 

9 

May 

101 

108 

115 

108 

107 

93 

44 

Port  Said,  Egypt 

45 

May 

96 

99 

102 

102 

99 

93 

45 

Tashkent,  USSR 

19 

May 

93 

100 

105 

101 

91 

74 

46 

Tefia,  Canary  Is, 

6 

June 

95 

99 

100 

95 

91 

84 

Bs  Climates 

47 

Ankara,  Turkey 

26 

May 

'  90 

94 

98 

99 

90 

78 

48 

Baghdad, Iraq 

15 

May 

107 

112 

116 

115 

no 

97 

49 

Berbera,  Br.  Som, 

30 

May 

110 

117 

115 

114 

113 

105 

50 

Dakar,  Senegal 

16 

Jun 

104 

104 

103 

105 

103 

97 
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TABLE  IV  (CONT.) 
AVERAGE  INTERNAL  HIGH  TEMPERATURES  - 

6  HOTTEST  MONTHS 

Sta 

No 

Name 

Yrs. 

Data 

1st 

Mo 

1 

2 

3 

**4 

5 

6 

51 

Miles  City,  Mon. 

39 

Apr 

76 

86 

95 

102 

100 

84 

52 

Kabul,  Afgh. 

8 

May 

95 

101 

105 

103 

98 

84 

53 

Jerusalem,  Israel 

18 

May 

91 

97 

101 

98 

96 

90 

54 

Haifa,  Israel 

16 

May 

97 

100 

102 

103 

99 

93 

55 

Flagstaff,  Ariz. 

20 

May 

86 

95 

98 

96 

89 

77 

56 

San  Diego,  Cal. 

59 

May 

85 

88 

91 

91 

88 

81 

57 

Santa  Fe,  N.M, 

54 

May 

86 

94 

97 

95 

89 

75 

58 

Ulan  Bator,  Mong. 

12 

May 

76 

88 

95 

86 

73 

54 

Ca  Climates 

59 

Brisbane,  Aus. 

52 

Oct 

94 

96 

99 

98 

95 

96 

60 

Buenos  Aires,  Arg. 

23 

Nov 

91 

98 

100 

98 

92 

82 

61 

Chungking,  China 

13 

May 

95 

101 

108 

107 

97 

85 

62 

Durban,  S. Africa 

15 

Nov 

90 

92 

93 

93 

90 

90 

63 

Hankow,  China 

13 

May 

96 

103 

108 

108 

97 

87 

64 

Hongkong,  China 

50 

May 

92 

94 

96 

96 

93 

86 

65 

Jacksonville,  Fla. 

59 

May 

99 

103 

106 

105 

100 

90 

66 

Louisville,  Ky. 

58 

May 

91 

99 

103 

100 

93 

78 

67 

Luang  Prabang.Mala  12 

Mar 

99 

102 

102 

100 

98 

98 

68 

Milano,  Italy 

10 

May 

89 

96 

98 

96 

92 

68 

69 

Mobile,  Ala. 

30 

May 

98 

103 

106 

105 

99 

89 

70 

Nagasaki,  Japan 

35 

May 

91 

95 

102 

104 

95 

83 

71 

Norfolk,  Va. 

60 

May 

91 

97 

101 

98 

93 

80 

72 

Oklahoma  City,Okla.  9 

May 

95 

102 

107 

107 

98 

83 

73 

Osaka,  Japan 

65 

May 

90 

97 

102 

105 

96 

83 

74 

Shanghai,  China 

38 

May 

95 

99 

106 

106 

96 

86 

75 

Tokyo,  Japan 

35 

May 

89 

93 

100 

102 

92 

78 
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TABLE  IV  (CONT.) 

AVERAGE  INTERNAL  HIGH  TEMPERATURES  -  6  HOTTEST  MONTHS 


Sta 

No 

Name 

Yrs. 

^Data^ 

1st 

Mo 

1 

L  ~T 

T] 

4 _ 

5 

6 

76 

Venice,  Italy 

10 

May 

89 

95 

97 

96 

90 

73 

77 

Washington,  D.C. 

60 

May 

89 

95 

97 

96 

90 

78 

Cs  Climates 

78 

Algiers,  Algeria 

25 

May 

88 

96 

99  . 

99 

93 

82 

79 

Athenai,  Greece 

71 

May 

92 

98 

102 

100 

93 

81 

80 

Beirut,  Lebanon 

62 

May 

93 

97 

101 

102 

98 

88 

81 

Bizerte,  Tunisia 

28 

May 

89 

97 

102 

104 

96 

86 

82 

Capetown,  S. Africa 

19 

Nov 

86 

89 

90 

91 

88 

84 

83 

Erivan,  USSR 

17 

May 

92 

101 

106 

105 

93 

78 

84 

Famagusta,  Cyprus 

30 

May 

96 

102 

107 

108 

101 

92 

85 

Iraklion,  Crete 

21 

May 

90 

95 

98 

98 

93 

85 

86 

Istanbul,  Turkey 

18 

May 

85 

94 

97 

96 

88 

73 

87 

Lisbon,  Portugal 

75 

May 

85 

92 

94 

94 

88 

76 

88 

Madrid,  Spain 

30 

May 

88 

96 

102 

99 

91 

75 

89 

Napoli,  Italy 

10 

May 

90 

97 

100 

99 

92 

78 

90 

Perth,  Australia 

43 

Nov 

91 

95 

98 

97 

93 

84 

91 

Nicosia,  Cyprus 

39 

May 

96 

102 

107 

107 

101 

89 

92 

San  Francisco,  Cal, 

40 

May 

82 

86 

85 

83 

83 

78 

93 

Tangier,  Tan. 

35 

May 

87 

92 

96 

95 

88 

82 

94 

Valencia,  Spain 

26 

May 

89 

93 

98 

97 

93 

81 

95 

Valletta,  Malta 

90 

May 

88 

95 

99 

100 

94 

83 

96 

Valparaiso,  Chile 

27 

Nov 

87 

88 

98 

97 

92 

78 
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TABLE  IV  (CONT.) 

AVERAGE  INTERNAL  HIGH  TEMPERATURES  -  6  HOTTEST  MONTHS 


Da  Climate « 


97 

Bucaresti,  Romania 

41 

May 

90 

97 

101 

100 

88 

73 

98 

Cagliari,  Sardinia 

10 

May 

88 

96 

101 

99 

93 

82 

99 

Pusan,  Korea 

36 

May 

86 

92 

98 

100 

91 

80 

100 

Roma,  Italy 

10 

May 

90 

98 

102 

101 

94 

80 

101 

Toulon,  France 

49 

May 

87 

94 

95 

95 

87 

76 

H  Climates 

102 

Addis  Ababa,  Eth, 

12 

Jan 

90 

92 

94 

96 

95 

90 

103 

Mexico  City,  Mexico 

7 

Apr 

95 

96 

96 

94 

94 

93 

Cb  Climates 

104 

Bordeaux, France 

41 

May 

86 

93 

94 

95 

87 

74 

105 

Cherbourg,  France 

16 

May 

79 

82 

85 

85 

77 

62 

106 

Frankfurt,  Germany 

50 

May 

84 

89 

92 

89 

79 

60 

107 

Greenwich,  England 

30 

May 

80 

86 

89 

86 

77 

58 

108 

Koln,  Germany 

50 

Apr 

74 

83 

87 

89 

86 

77 

109 

Paris,  France 

24 

May 

84 

89 

91 

91 

86 

65 

110 

Renfrew,  Scotland 

30 

Apr 

70 

75 

82 

85 

79 

71 

111 

Rotterdam,  Neth, 

30 

May 

80 

84 

87 

83 

75 

60 

112 

Wellington,  N,Z, 

65 

Nov 

82 

85 

86 

85 

80 

70 

Db  Climates 

113 

Kiev,  USSR 

10 

May 

86 

93 

94 

90 

78 

58 

114 

Nemuro,  Japan 

65 

May 

73 

78 

85 

87 

80 

66 

115 

Vladivostok,  USSR 

14 

May 

76 

83 

89 

91 

80 

64 

Dc  Climates 

116 

Archangel,  USSR 

23  ' 

Apr 

37 

66 

75 

80 

79 

55 

117 

Churchill,  Manitoba 

15 

May 

63 

71 

84 

76 

59 

35 

118 

Fairbanks,  Alaska 

13 

Apr 

58 

75 

86 

88 

79 

59 

119 

Goose  Bay,  Newf. 

6 

May 

70 

79 

87 

81 

64 

48 

120 

Petropavlovsk,USSR 

7 

May 

67 

73 

77 

78 

68 

49 
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COMPUTING  WATER  GAIN  (GENERAL) 


Let 

Ti 


Pa 

Pv 

Wa 

Ww 

R 


=  temperature  at  which  water  stops  being  inhaled  on 
any  given  day. 

=  temperature  at  which  water  begins  to  be  inhaled. 

=  pressure  at  which  water  is  inhaled. 

=  partial  pressure  of  air. 

=  partial  pressure  of  water  vapor. 

=  weight  of  dry  air  inhaled  per  cubic  foot  of  original 
volume  on  any  given  day. 

=  weight  of  water  inhaled  per  cubic  foot  of  original 
volume  on  any  given  day. 

=  gas  constant  for  dry  air,  53,  35  ft-lb/lb/°R. 

=  specific  humidity  of  air  inhaled. 


Then,  simple  manipulation  of  the  basic  P-V-T  relationship  for 
a  constant  voluihe  and  constant  inhale  pressure,  gives 


Hl 

R 


(27) 


and 


s  W0 


(28) 


*  This  assumes  that  the  oil  canning  effect  is  negligible  compared  to 
total  volume. 
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Now,  by  definition  ( 


p„  p,  p,  h  r 

8'  1.61  Po‘  1.61  (P,  -p  )“  t.SI  F»  (|_4J1)-  '«'P, 

P, 

Since  the  daily  amount  is  small,  it  is  desirable  to  express 
Ingrams.  Usual  units  for  pg  and  Pj  are  inches  of  mercury 
that  F  is  a  dimensionless  number.  Combining  (27)  and  (28) 
have  Ww  »  in  grams,  for  a  given  day,  by 


<29/ 

Ww 

such 

we 


Ww  =  C.ocr  (30) 

where  Cj  =  758.32588. 

For  n  days,  total  weight  of  water  inhaled  is,  obviously,  given  by 


w„, .  c,  |  <,rt  (51» 

where  the  subscripts  i  denote  the  days  on  which  water  was  actually 
inhaled.  On  free  breathing  containers  this  is  every  day  but  this  does 
not  necessarily  hold  when  check  valves  prevent  free  passage  of  air,  a 
fact  which  will  become  important  later. 
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FREE  BREATHING  CONTAINERS 
(GENERAL  SURVEY) 


Average  Water  Gains 

In  free  breathing  container*,  Tj  °  Tj  and  Tj  *  the  ambient 
low  reached  at  night.  In  computing  I*  .  the  temperature  and  relative 
humidity  are  constantly  changing.  Composition  of  the  inhaled  air  may, 
however,  be  characterised  by  taking  the  average  of  the  high  and  low 
relative  humidities  and  pg  at  the  average  of  Tj  and  T^. 

As  noted  previously,  tj  was  determined  by  hand  for  each  mid¬ 
month  of  the  six  hottest  months  for  120  stations.  This  tj  was  taken  as 
the  average  for  the  month  and  the  indicated  computations  were  made  for 
each  of  the  120  stations. 

Detailed  computing  results  for  one  station  are  shown  in  table  V . 
Total  water  gains  for  all  stations  are  shown  in  figure  3. 

It  is  particularly  interesting  to  note  that  the  monsoon  climate 
(Am)  appears  most  severe  even  though  the  daily  temperature  change  is 
small.  Tropical  rainforest  (Af)  and  tropical  savanna  (Aw)  are  a  close 
second  and  third.  Further,  note  the  high  position  occupied  by  deserts 
(Bw).  In  fact,  although  humid  subtropical  (Ca)  —  typically  southeastern 
United  States  --is  shown  in  fourth  position  because  of  a  higher  mean, 
the  extreme  high  values  for  Bw  are  comparable  to  the  more  severe 
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AVERAGE  WATER  GAIN  PER  CUBIC  FOOT  AT  OKLAHOMA  CITY 
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DESICCANT  BREATHER  SERVICE  LIFE 


The  probable  service  life  of  a  desiccant  breather  depends  on 
the  water  flow  in  ground  storage,  the  number  of  flights,  the  extent  of 
reactivation  by  dry  container  air  flowing  out,  the  amount  of  desiccant, 
and  the  water  sorptive  capacity  of  the  desiccant.  The  latter,  in  turn, 
is  a  function  of  the  point  at  which  the  desiccant  is  considered  exhausted. 

A  general  expression  for  service  life  in  months  per  gram  of 
desiccant  per  cubic  foot  of  container  volume  (Mj)  is  easily  derived  as: 


M. 


M,  W, 

(I  -C)(W„,  +nW,) 


where 


(32) 


M2 

Wwt 

W2 

c 

n 

W3 


=  months  of  exposure.  In  this  study  M2  =  6. 

=  total  water  gain  in  time  period  M2 

=  percentage  water  storage  capacity  of  desiccant  at 
its  end  point,  expressed  decimally. 

a  desiccant  reactivation  factor. 

=  number  of  airplane  descents. 

=  grams  of  water  inhaled  per  descent  per  cubic  foot 
of  container. 


In  preliminary  conversations  leading  up  to  this  work  (21)^  it  was 
agreed  that  the  logical  endpoint  should  be  when  influent  air  reaches  40 
percent  R.  H.  Thus,  using  Type  I,  Grade  H,  desiccant  conforming  to 
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specification  (22)^  =  0.215.  Reworking  the  Peterman  and  Nelson 

data  given  in  Appendix  I  for  a  free  breathing  container  gives  W3  = 
0.426  grams  H2O  per  cubic  foot  of  container  volume. 

Gelber'  1  ha8  made  preliminary  measurements  of  reactivation 
values.  Where  exhaled  air  was  at  10  percent  RH,  C  =  0.312.  At  0 
percent  RH,  C  =  0.  527.  It  is  reasonable  to  assume  that  a  drawdown 


VALVED  CONTAINERS 


The  General  Problem 


Ty  =  temperature  at  which  the  vacuum  relief  valve  opens. 

Tp  =  temperature  at  which  the  pressure  relief  valve  opens. 

Pp  =  pressure  at  which  vacuum  relief  valve  opens. 

=  pressure  at  which  pressure  relief  valve  opens. 

^low  =  ambient  low  temperature. 

Ti  and  T2  are  already  defined  in  equation  (27) . 


Then,  assuming  constant  volume,  as  before,  we  have 


t,  =  -jr-V  T. 

rv 


If,  and  only  if,  Tjow  <  Ty  on  a  given  day,  the  container  will 
gulp  water  on  that  day  and  equations  (27)  and  (31)  will  come  into  play  with 
T2  =  T^ow.  Note,  howeVer,  that  when  Tjow  <  Ty,  the  gas  law 
indicates  a  new  Tp  for  the  next  day  by 

T  ■  -£■  T*.  (35) 

p 

If,  however,  on  this  next  day,  T^  >  Tp,  then  we  have  a  new 
Tv  from 


Tv  ft  Ti 


Thus,  it  is  clear  that  behavior  of  a  valved  container,  on  any 
given  day,  is  dependent  on  the  past  history  of  the  container  in  a  specific 
climate.  It  is  evident  also  that  reliance  on  monthly  average  tempera¬ 
tures  to  determine  water  intake  over  a  period  of  time  simply  will  not 
give  a  true  picture  of  the  probabilities  and  will,  therefore,  fail  to  pro- 
vide  the  information  sought  herein. 

On  valved  containers,  we  are  interested  in  the  outcomes  (water 
gains)  of  a  series  of  possibly  highly  variable  events.  In  order  to  arrive 
at  usable  design  values,  we  must  achieve  a  statistical  distribution  of 
these  possible  outcomes.  The  next  section  will  discuss  the  statistical 
problem  in  some  detail  and  demonstrate  how  a  usable  answer  can  be 
achieved. 

Before  turning  to  the  statistical  problem,  it  is  worthwhile  to 
define  the  valving  arrangements  actually  considered  out  of  the  multitude 
of  those  possible.  The  following  configurations  were  considered  since 
they  are  representative  of  feasible  commercially  available  combinations: 


*  It  should  be  noted,  further,  that  even  on  free  breathing  containers, 
we  have  not  yet  found  the  distribution  about  the  mean  values  com¬ 
puted. 
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Relative 

Relative 

Pressure 

Vacuum 

Configuration 

Relief,  psig 

Relief,  psig 

A 

None 

None 

B 

+  2 

-  1 

C 

+  1 

-  1 

D 

+  1 

-  1/2 

E 

+  1/2 

-  1/2 

Note  that 

Py  =  Pj  in  equations  (27)  and 

(29)  and  that  this 

pressure  is  dependent  upon  the  station  altitude.  Correction  for  alti- 

tude  is  readily  obtained  from 

p 

'  p« 

P,' 

(37) 

where 

=  sea  level  pressure  at  which  water  is 

inhaled. 

The  ratios  Pv/Pp  and 

P  /P  are  constant  for  each  valve  con 
p  v 

figuration.  Values  of  p{  and  these  ratios  are  tabulated: 

Valve 

Configuration 

p/ 

(in,  Hg.) 

pv/pp 

P  /P 

p  V 

A 

29.92 

B 

27.  88 

. 82024 

1.21915 

C 

27.88 

.  87234 

1. 14634 

D 

28.90 

.90426 

1. 10588 

E 

28.  90 

.93369 

1.07024 
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The  Statistical  Problem 


As  might  be  expected  from  the  dependence  of  I^av  on  a  sine 
function  —  equations  (24),  (25)  and  (26)  --  mean  high  and  low  temper¬ 
atures  can  usually  be  matched  to  a  sine  function.  Extensive  study  re¬ 
veals  that  the  full  period,  2  ir  ,  covers  the  entire  year  in  the  temper¬ 
ate  zones.  At  the  equator,  on  the  other  hand,  the  full  year  encompasses 
$ 

two  periods. 

The  temperature  data  in  the  world-wide  tables  U&)  consist  of 

% 

(among  other  data),  the  mean  high  and  low  temperatures  for  each  month 
and  the  number  of  years  the  data  have  been  collected.  A  typical  sine 
match  for  one  city  is  shown  in  figure  5.  The  matching  technique  makes 
the  mean  of  the  fitted  curve  equal  to  the  mean  of  the  empirical  data  and 
the  area  under  the  fitted  curve  equal  to  the  area  under  the  data  curve. 
Note  that  there  is  a  one-month  phase  angle  lag  behind  the  insolation 
values  derived  from  equations  (20),  (25)  and  (26)  which  is  in  accord  with 
theoretical  expectancy.  ** 


*  As  location  approaches  the  tropics  of  Capricorn  or  Cancer,  the 
second  period  reduces  in  both  amplitude  and  time  duration.  In 
theory  this  should  be  taken  into  account  in  any  computing  program. 
In  practice,  here,  it  was  not  and  errors  were  minimized  by  taking 
stations  as  far  north  as  practicable.  Since  fluctuations  are  small 
in  tropical  climates  anyhow,  absolute  error  is  small. 

##  Once  again,  it  is  pertinent  to  note  that  heat  gain  by  the  atmosphere 
is  primarily  the  result  of  reradiation  by  the  earth  which  is  a  for¬ 
midable  heat  sink. 


MONTHS 


Figure  5  -  Temperature  at  Mobile,  Ala. 
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Wcsi 


Once  a  satisfactory  sine  curve  match  is  achieved  from  the 
tabulated  data,  then  the  probable  mean  high  and  low  temperature  for 
the  *  ith  day  and  the  mean  of  the  extreme  high  and  low  temperatures 
for  that  day  are  available  from 

Xj  *  x  ±  A  sin  cdj  (38) 

where 

x.  =  expected  mean  value  for  the  i  th  day 

x  =  annual  mean  value 

A  =  amplitude  of  the  sine  curve,  i.  e. ,  maximum  value 

of  x  above  or  below  x 

c  =  correction  factor  defined  in  equation  (26) 

d^  =  the  number  of  days  after  the  day  on  which  x  occurs, 

i.  e.  ,  after  April  22. 

Thun  we  have  a  pair  of  values  representing  the  mean  daily  high 
or  low  temperature  on  any  given  day.  We  have,  further,  the  number  of 
years  of  data  for  each  station  from  the  world-wide  tables  which 
defines  the  sample  size,  n.  Now,  we  have  it  on  good  authority  (^4) 
the  actual  temperature  is  normally  distributed  about  the  mean  value 
after  removal  of  trend  values.  Thus  we  are  able  to  use  extremal  sta¬ 
tistical  theory  to  determine  the  number  of  standard  deviations,  n(s), 
represented  by  the  mean  and  extreme  value  pair  previously  derived. 
Figure  6,  taken  from  Graph  4.  2.  2.  (4)  of  Gumbel  (25)^  t^e  means  for 
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n(s)  vs  Sample  Size  for  Mean  Extreme  Values 
(after  Gumbel) 


accomplishing  this  statistical  operation.  Note  that  n(s)  is  a  constant 
for  each  station  analyzed. 

Now,  let 

xi  =  mean  high  (low)  temperature  for  the  year. 

X2  =  mean  of  the  extreme  high  (low)  temperatures  for 
the  year. 

Ai  =  amplitude  of  mean  high  (low)  temperature  curve 
for  the  year. 

A2  =  amplitude  of  the  average  extreme  high  (low) 
temperatures  for  the  year. 

xu  =  mean  high  (low)  temperature  on  the  i  th  day 

x2i  =  averaf5*  extreme  high  (low)  temperature  on  the  ith 

day*. 

Then,  the  size  of  the  standard  deviation  in  degrees,  sj,  for  the 
high  or  low  temperature  on  the  i  th  day  is  given  by 

(39) 

Thus,  on  the  i  th  day,  equation  (38)  gives  us  the  mean  high 
(low)  for  the  day  and  equations  (39)  gives  the  dimensioned  standard  de¬ 
viation  for  that  day.  These  values,  coupled  with  I^y  *or  each  day 
determined  from  equation  (24)  constitute  a  mathematical  model  of  the 
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day-to-day  internal  temperature  variation  inside  a  container. 

It  is  now  possible  to  construct,  for  any  given  day,  a  reasonable 
high  (low)  temperature  from  the  relation 


t 


LOW 


=  X;  +  CT  Sj 


(40) 


where 


<7  =  a  signed  random  normal  deviate  taken  from  a  table  of 

random  normal  deviates  based  on  a  truly  random  col¬ 
lection  of  numbers,  e.  g. ,  the  Rand  Table 


With  tQ  thus  determined,  it  is  possible  to  compute  t^  and 
compare  the  valve  settings.  With  tjQW  it  is  possible  to  determine 
whether  or  not  the  valve  opens  and,  if  so,  how  much  water  is  inhaled. 
Repeating  this  process  365  times  we  have  a  year's  water  gain  for  each 
valve  configuration.  Repeat  for  n  years  at  each  station  and  we  will 
ha,'.-!  a  statistical  distribution  of  water  gains  covering  a  sample  of  size 
n.  From  this  distribution  it  is  then  possible  to  draw  valid  statistical 
inferences  concerning  the  water  gain  to  be  expected  in  specific  climates. 

It  should  be  particularly  noted  that,  although  high  and  low  temper¬ 
atures  are  randomly  distributed  about  their  means  in  a  large  sample, 
the  high  and  low  temperature  on  a  given  date  are  not  completely  inde¬ 
pendent  events.  Similarly,  there  is  usually  about  a  three  day  persistence 


*  This,  of  course,  is  using  the  Monte  Carlo  method  in  one  of  its  simpler 
forms.  Justifying  the  validity  of  using  the  results  of  this  form  of  simu 
lating  reality  is  outside  the  scope  of  this  report. 
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effect  for  any  given  weather  cycle.  In  order  to  simulate  these  practical 

> 

constraints  on  the  chances  of  extreme  variation,  use  is  made  of  the  cor¬ 
relation  formula  (27) 


o r  =  p  <r 1  +  <r  Vl  ~  P2 


(41) 


wiiere 

=  correlated  new  deviate 
=  previous  deviate 

=  new  random  deviate  taken  from  the  tables 
=  degree  of  correlation 

Day  and  night  correlation  is  set  at  90  percent.  Correlation  be¬ 
tween  successive  days  is  set  at  75  percent  for  three  days,  after  which  a 
completely  new  day  random  deviate  is  taken. 

Equations  (1)  through  (41)  thus  define  a  mathematical  model  of 
the  effects  of  a  climate  at  any  selected  location.  While  each  computation 
is  relatively  simple,  there  are  a  large  number  of  them.  Obviously  use 
of  a  computer  is  necessary. 

A  general  digital  computer  arrangement  is  shown  in  figure  7. 

The  climate  conversion  and  storage  processes  are  shown  in  figure 
8  while  the  operation  of  generating  a  daily  high  internal  and  low  ambient 
temperature  is  shown  in  figure  9.  The  unit  which  decides  whether  or  not 
water  is  gulped  on  a  given  day  is  diagrammed  in  figure  10.  Water  gain 
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(T)  Perform365  times 

Perform  n  times  and  stop 


Figure  7 


Basic  Computer  Routine 
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Generator 


Water  Gain  Calculator  and  Print  Out 


f 

computation  and  print  out  is  shown  in  figure  11. 

Actual  programming  and  detailed  computations  were  undertaken 
by  the  Space  and  Information  Systems  Division  of  North  American  Avi¬ 
ation,  under  the  direction  of  Messrs.  Masaji  Hatae  and  Elliott  Kleinman. 

Since  tj  is  dependent  on  tQ,  the  question  will  naturally  arise  con¬ 
cerning  why  not  simply  operate  statistically  on  tj  directly.  The  Under¬ 
lying  reason  for  not  doing  so  is  twofold: 

a)  Transforming  tQ  to  t^  involves  a  non-linear  relation.  The 

>  basic  operators  are  ^  and  j>t  which  vary  as  shown  in 

figure  1,  producing  a  throttling  effect  on  maximum  and  minimum 
tj.  That  is  to  say,  the  variation  in  tj  is  less  than  the  variation 
in  tQ. 

b)  The  only  reasonably  statistically  valid  estimate  attainable  is 
the  daily  mean  and  standard  deviation  of  tQ. 

Selecting  the  Stations  to  be  Analyzed 

While  the  individual  computations  diagrammed  in  figure  7  through 
11  are  each  relatively  simple,  there  are  a  very  large  number  of  them. 

For  obvious  reasons,  it  is  necessary  to  limit  the  number  of  stations. 

Detailed  input  data  on  the  stations  analyzed  are  shown  in  table  VI. 
The  following  reasoning  governed  selection: 

a.  Since  we  are  primarily  interested  in  being  able  to  ship  and 
store  anywhere,  low  water  gain  climates  can  be  ignored.  Based  on  the 


TABLE  VI 


Statistical  Data  on  Stations  Analyzed  by  Computer 


information  in  figure  3,  H,  Cb,  Db  and  Dc  climates  were  dropped. 

b.  The  Da  climate  —  humid  continental,  warm  summer  — 
appears  adequately  bracketed  by  Ca  and  Cs  climates  so  that  stations 
from  the  latter  two  cover  the  design  problem. 

c.  Stations  are  taken  from  the  higher  valued  gains  in  each 
climate  group  because  we  are  seeking  a  solution  to  the  design  problem. 

d.  Temperature  sine  curves  at  a  given  station  should  match  the 
empirical  data  acceptably.  The  criterion  of  acceptability  was  P=0.  99 
for  mean  curves  and  P=.  95  for  extreme  curves  using  the  Chi-squared 
statistic.  With  ten  degrees  of  freedom,  anticipated  values  of  Chi- 
squared  for  these  probabilities  are  2.558  and  3.940  respectively (28)., 
All  curves  selected  exceed  99  percent  significance  except  the  extreme 
value  curves  for  Famagusta,  Cyprus  which  exceed  95  percent  signifi- 
cance. 


#  Calcutta  showed  a  Chi-squared  value  of  5.  703  on  mean  extreme 
highs  while  Rangoon  showed  6.  375  and  6.  6111  on  mean  high  and 
mean  extreme  high.  Here,  ambient  temperature  drops  at  a  time 
when  theoretically,  it  should  be  rising,  apparently  because  of  the 
onset  of  the  rainy  season.  Because  of  poor  match,  these  stations 
were  not  analyzed  by  computer. 
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RESULTS  OF  COMPUTER  RUNS 


Unvalved  Containers 

Results  of  computer  runs  for  unvalved  containers  are  shown 
in  tables  VII  and  VII.  These  tables  show,  for  six  months  and  one  year 
respectively,  the  mean  gain  (X),  the  sample  standard  deviation  (S)  and 
the  various  constants  needed  to  compute  a  distribution  curve  based  on 
the  sample. 

Since  the  a  priori  distribution  function  was  not  known,  the 


Pearson  system  --  summarized  in  the  literature  (29*30)^  Q£ 
sification  was  used.  In  this  system  the  following  properties,  based  on 
the  first  four  moments,  of  the  empirical  distribution  were  obtained: 


a.  The  mean,  defined  by 

X-jrZXj 


(42) 


b.  The  standard  deviation 


S« 


/ 


S*i  -2 

- X 


(43) 


c.  The  skewness 


0  3  * 


2(X|-x)S 


ns  3 


(44) 


d.  The  peakedness 

z(xi'x)4 


ns 


(45) 
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Annual  Water  Gains 
Free  Breathing  Container 
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In  the  Pearson  system  of  classification,  the  first  step  is  to 
compute  the  reduced  critical  function, 

F  a  2a  -  3af  -  6  (46) 

I  4  3 

As  shown  in  the  tables  the  computed  value  of  this  function  is 
negative  in  every  case  but  one  and  a^  is  positive  in  every  case. 
Thus  each  distribution  curve  except  that  for  San  Juan  is  a  Pearson 
Type  I,  defined,  with  the  origin  at  the  mode,  by 


i,  8  j  [i  -Dr)  (52) 

i2  =  i  -  i,  (53) 

J I 

m|S  -i  (  R  -  2)  (54) 

m  j  +•  m g  =  R  -  2  (55) 


The  Type  I  distribution  function  is  limited  in  both  directions 
and  is  normally  ske'ved.  If  mj  and  are  approximately  equal 

the  curve  is  nearly  symmetrical.  If  mj  and  m^  are  not  small, 
it  tails  off  at  both  ends  but  if  both  are  small  it  rises  abruptly  at  both 
ends.  If  mj  is  negative  the  curve  is  J  -  shaped  while  if  both  mj  and 
are  negative  the  curve  is  U-shaped,  lj  and  1^  are  the  distances  on 
either  side  of  the  mode  defining  the  limits  of  the  curve. 

In  the  semi-yearly  curves,  no  data  point  exceeds  I2  but  this 
is  not  true  for  the  yearly  values,  thus  indicating  a  poor  match.  Both 
tables  VII  and  VIII  show  the  maximum  value  recorded  and  the  number  of 
standard  deviations  (t)  this  value  represents.  It  is  apparent  that  maxi¬ 
mum  t  is  very  high  --  in  fact  these  values  often  could  be  rejected  under 
Chauvenet's  criterion.  They  have  not  been  rejected,  however,  because 
we  have  no  physical  basis  for  doing  so,  even  though  we  are  not  able  to 
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draw  valid  statistical  inferences  from  the  data. 


In  seeking  a  cause  for  the  difficulty  on  the  one  year  data,  the 
following  appears  reasonable: 

a.  Each  year  uses  exactly  the  same  sequence  of  random 
normal  deviates  at  each  station. 

b.  In  the  programming,  only  enough  random  normal  devi¬ 
ates  were  placed  on  cards  to  run  for  18  years.  Between  runs,  the 
cards  were  thoroughly  shuffled  to  produce  a  new  set  of  deviates.  In 
effect,  we  have  two  separate  samples  drawn  from  the  same  universe. 
For  the  analysis  in  tables  VII  and  VIII  the  computed  moments  were 
combined  in  conventional  fashion. 

c.  In  the  second  half  of  year  six  of  the  second  run  an  unu¬ 
sual  combination  of  random  deviates  must  have  occurred  since  Xm 
invariably  is  found  here.  The  resulting  high  water  gain  changed 
curves  which  were  normally  negatively  skewed  to  positive  skew  and 
also  increased  peakedness  significantly.  Generation  of  this  unusually 
high  number  --  which  represents  a  possible  climate  --  rendered  the 
data  non-homogeneous . 

d.  The  cure  for  non-homogeneous  data  from  a  presumably 
homogeneous  universe  is  only  found  by  taking  a  very  much  larger 
sample,  a  proceeding  not  possible  at  this  time. 

e.  In  effect,  we  gambled  --  from  early  check  runs  indicating 
stability  and  a  very  small  standard  deviation  --  that  a  small  sample 
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would  be  ample  and  lost.  Fortunately,  the  sample  was  sufficiently 
large  that  a  false  sense  of  security  could  be  avoided. 

It  is  interesting  to  note  that  the  single  wild  data  point  intro- 
duces  some  provocative  thoughts.  Specifically,  although  the  mean 
gain  for  the  tropical  and  monsoon  stations  is  the  highest,  in  accord¬ 
ance  with  expectancy,  the  highest  single  one  year  gain  occurs  at 
Las  Vegas,  a  desert  station.  Here  we  have  a  single  year  gain  of 
31.  15  grams  per  cubic  foot  whereas  maximum  gain  at  Bangkok  is 
22.25  grams  pe.r  cubic  foot.  This  reinforces  the  suspicion  that 
desert  climates  can  be  as  severe,  or  more  so,  from  a  water  gain 
standpoint  as  any  in  the  world. 

We  are  permitted,  however,  to  assume  that  the  mean  water 
gain  is  normally  distributed  about  the  true  mean  gain  and  the  standard 
deviation  of  this  distribution  is  a  function  of  the  sample  standard  devi¬ 
ations.  Thus,  for  Las  Vegas,  we  have  in  the  one  year  runst 
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Since  y  is  less  than  twice  s  (y),  the  differences  between 
X}  and  X^  are  probably  not  significant.  Hence,  we  can  use  the 
computed  mean  of  the  combined  sample  (7.  17)  with  good  confidence 
that  its  probable  error  does  not  exceed  +  .  318  grams  of  water  per 
cubic  foot  of  container  volume.  The  example  chosen  is  the  worst 
case. 

It  is  useful  to  speculate  briefly,  on  what  kind  of  distribution 
function  could  reasonably  be  expected  from  a  sufficiently  large  sample 
(say  n  200).  The  distribution  is  definitely  limited  to  the  left  (it 
cannot  be  less  than  zero)  and  may  or  may  not  be  limited  to  the  right. 
Positive  skewness  is,  apparently,  probable.  Accordingly,  a  Type  I 
function,  possibly  degenerating  into  a  Pearson  Type  III,  is  indicated. 
The  least  likely  distribution  is  the  handy  normal  distribution  function. 

The  single  exception  to  the  Type  I  functions  is  the  yearly  gain 
at  San  Juan  where  Fj  is  slightly  positive.  It  is  convenient  here,  to 
assume  a  Type  III  distribution  function,  defined,  with  the  origin  at  the 
mode,  by 
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(58) 


(59) 


Curve  is  limited  on  the  left  if  a-j  is  positive.  The  curve  is  usually 
bell  shaped  but  becomes  J  shaped  when  is  greater  than  4. 

Table  VIII  lists,  therefore,  only  lj. 

It  is  interesting  to  compare  the  hand  computed  six  months 
values  with  the  machine  computed  values.  The  results  are  shown  in 
table  IX. 

A  detailed  investigation  of  the  reasons  for  the  large  differences 
is  not  necessary.  Suffice  it  to  say  that  the  machine  computing  method 
gives  greater  confidence  and,  for  any  specific  problem,  is  the  preferred 
method. 

Taking  the  mean  values  from  table  VII  and  substituting  into 
equation  (33),  we  obtain  the  mean  storage  lives  per  gram  of  desiccant 
in  the  breather  shown  in  table  X.  It  is'apprent  that  wide  fluctuation  in 
lives  is  possible  depending  upon  the  number  of  flights  encountered  and 
the  location.  Recasting  table  X  in  terms  of  English  units,  we  have  table 
XI.  It  is  a  simple  slide  rule  operation  to  convert  the  values  of  table  XI 
to  any  other  loadings  which  may  be  desired. 
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Table  EX 


Comparison  of  Hand  and  Monte  Carlo  Results 


Station 

Hand 

Machine 

Madras 

6.  12 

10.64 

Manila 

6.  55 

10.31 

Berbera 

5.  58 

8.63 

Famagusta 

4.  18 

6.44 

Bahrein 

5.82 

10.72 

Bangkok 

6.58 

11.78 

Washington 

5.38 

'  5.57 

Mobile 

4.99 

8.46 

Las  Vegas 

2.32 

2.40 

< 

Belize 

6.48 

10.83 

Miles  City 

5.55 

3.93 

San  Juan 

4.82 

7.21 

• 

Percentage 

Diff 

+ 

75 

+ 

58 

+ 

54 

+ 

55 

+ 

84 

+ 

81 

Months  of  Life  Per 


Ounce  of  Desiccant  Per  Cubic  Foot 


Based  upon  perusal  of  table  XI,  it  is  recommended  that  the 
standard  breather  loading  for  world  wide  distribution  be  set  at  one 
ounce  per  cubic  foot  of  free  air  volume. 

A  safety  factor  working  for  better  service  lives  than  those 
indicated  is  the  desiccating  ability  of  the  breather  tubes  themselves. 
While  no  complete  study  of  this  effect  has  been  made  —  and  in  the 
author's  opinion  one  should  --  Zerr  (31)  found  that  five  1/4  x  2-1/2 
inch  tubes  in  parallel  would  reduce  relative  humidity  of  80°F.  air 
from  90  to  76  percent  at  flow  rates  as  high  as  0.4  cubic  feet  per 
minute. 

At  the  same  time,  however,  the  computations  made  here 
assume  that  the  only  source  of  water  is  outside  the  container.  Water 
trapped  in  the  container  should  be  purged  dynamically,  if  possible. 
Otherwise  this  water,  even  if  initially  captured  by  a  static  draw  down 
charge,  will  eventually  be  shared  by  the  breather  desiccant  with  ap¬ 
propriately  shortened  life. 

Valved  Containers 

Turning,  now,  to  the  check  valves,  we  arrive  at  some  reason¬ 
ably  unexpected  results.  In  early  check  runs  it  was  found  that  annual 
water  gains  were  very  low.  Accordingly,  the  program  was  modified 
to  count  the  number  of  times  a  valve  opened  in  each  climate.  Where 
the  number  of  openings  was  small,  the  probability  of  a  given  number  of 
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openings  was  computed  by  the  Poisson  distribution. 


P(n) .  fils! 

where 

n  =  number  of  openings 

n  =  average  number  of  openings  in  the  sample 

Actually,  of  course,  we  are  dealing  with  a  binomial  problem 
(open  -  not  open).  If  the  number  of  trials  is  large,  and  the  number  of 
openings  is  fairly  high,  we  can  estimate  the  mean  number  of  openings 
and  the  standard  deviation  in  the  usual  way  and  determine  the  2<j 
number  of  openings  (representing  the  95%  probable  number  of  openings) 
using  the  normal  distribution  function.  These  two  methods  were  used 
in  the  detailed  analyses  for  each  valve. 

+2,  -1  valve 

This  valve  practically  never  opens.  It  opened  twice  in 
Washington,  D.  C.  in  36  years  and  at  no  other  location.  Using  formu** 
la  (60),  we  find  the  following  probabilities  for  Washington,  D.  C. 

P  one  opening  is  .  0523 

P  two  openings  is  .0016 

p  three  openings  is  33  x  10’^ 
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The  mean  of  the  two  water  gains  reported  is  0.  754  x  10"^ 
grams  per  cubic  foot,  a  negligible  quantity.  Thus  the  quantity  of  ad¬ 
ditional  desiccant  which  may  be  required  is  completely  controlled  by 
the  flight  syllabus  and  is,  from  Appendix  I,  .  084  units  per  cubic  foot 
per  flight,  i.e. ,  0,  504  grams  of  water  per  cubic  foot  per  flight. 
MIL-P-116  required  1.2  units  per  cubic  foot  (on  an  empty  container 
basis)  with  a  total  water  capacity  of  about  7.2  grams.  Assuming 
sealing  at  standard  conditions  we  will  have  about  .  0648  grams  of 
water  per  cubic  foot  in  the  air  (ignoring  sorption  on  the  container 
walls).  Thus,  it  would  require  about  14  flights  to  saturate  existing 
desiccant. 

On  the  basis  of  the  foregoing,  it  is  concluded  that  a  +2,  -1 
valve  needs  no  additonal  desiccant  beyond  the  MIL-P-116  formula. 

41,  -1  valve 

This  valve  opens  rarely.  Total  number  of  openings  found  at 
Las  Vegas  is  17  with  mean  of  .472  per  year  and  S  =  .  784/year. 
Probability  of  one  opening  is  .  312,  of  two  is  .  038.  Thus  normal  de¬ 
sign  value  is  two  openings  per  year.  Mean  gain  per  opening  is  0.478 
x  10  grams  per  cubic  foot,  for  a  total  design  gain,  with  96.2%  con¬ 
fidence,  of  .  856  x  10“^  grams  per  year,  again  a  negligible  quantity. 

On  this  basis,  additional  desiccant  for  this  valve  is  not  re¬ 


quired. 


+1,  -1/2  valve 


This  valve  breathes  in  Belize,  Washington  and  Las  Vegas. 

The  number  of  gulps  at  Belize  was  two  in  36  years.  Proba¬ 
bilities  are,  therefore,  the  same  as  previously  computed.  The  mean 
size  of  the  gulp  was  .2389  x  10“^  grams  per  cubic  foot,  still  a  negli¬ 
gible  amount. 

The  mean  number  of  gulps  per  year  at  Las  Vegas  is  3.  583 
with  standard  deviation  2.476  (maximum  number  was  11).  The  de¬ 
sign  number  of  gulps  is,  therefore,  nine  per  year  with  98.  57%  con¬ 
fidence.  The  mean  gain  per  gulp  is  .  6624  x  10"®  grams  per  cubic 
foot.  Thus,  the  design  annual  gain  is  .  5962  x  lO"^  grams  per  cubic 
foot.  In  the  light  of  the  reserve  capacity  required  by  MIL-P-116, 
this  gain  is  negligible. 

Washington  gulped  twice,  hence  probabilities  are  as  before. 
Mean  gulp  size  was  .3316  x  10”®  grams,  which  may  be  neglected. 

Hence,  it  may  be  concluded  that  this  valve  configuration  is 
usable  without  additional  desiccant 

+ 1/2,  -1/2  valve 

This  valve  breathes  in  all  cities  except  Manila,  San  Juan 
and  Bahrein.  Breathing  probabilities  are  summarized  in  table  XII. 

In  all  cases,  the  design  gain  is  negligible. 

This  last  statement  deserves  some  expansion.  The  design 
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gain  for  Las  Vegas,  where  breathing  is  by  far  the  most  prevalent, 
is  .  0663  grams  of  water  per  cubic  foot  per  year.  Since  MIL-P-116 
requires  1.2  units  of  desiccant  per  cubic  foot  and  since  each  unit 
will  hold  six  grams  of  water,  it  would  take  something  over  108  years 
for  this  infiltration  to  exhaust  this  desiccant  capacity  assuming  a 
completely  dry  container.  Even  if  we  assume  that  95  percent  of  the 
desiccant  capacity  is  consumed  on  initial  draw  down,  the  container 
life  is  5.  7  years,  which  is  far  beyond  most  current  military  storage 
objectives  and  is  certainly  in  excess  of  an  assumed  six  months  content 
reinspection  cycle. 

Appendix  I  does  not  compute  gains  for  this  valve  con¬ 
figuration.  It  is,  however,  somewhat  greater  than  the  free  breathing 
unit  and  somewhat  less  than  the  +2,  -1  valve.  For  the  latter,  we 
already  concluded  that  effect  of  five  flights  on  storage  life  was  negli¬ 
gible  and  so  the  same  conclusion  applies  here. 

Hence,  we  conclude  that  the  +1/2,  -1/2  valve  needs  no  ad¬ 
ditional  desiccant. 
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Table  XII 


Breathing  by  +1/2,  -1/2  Valves 


City 

1 

2 

3 

4 

5 

6 

7 

Madras 

2 

.056 

-- 

1 

.9477 

.  1691x10 

. 1691x10*2 

Manila 

0 

m  mm 

-- 

0 

-- 

-- 

0 

Berbera 

2 

.056 

mm  mm 

1 

.9477 

.8322xl0“3 

.8322xl0-3 

Famagusta 

11 

.306 

m  mm 

2 

.9657 

.2127x10-2 

.4254x10-2 

Bahrein 

0 

-- 

-- 

0 

— 

m  mm 

0 

Bangkok 

1 

.028 

1 

.9728 

,2366xl0"2 

.2366x10-2 

Washington 

4 

.  Ill 

-- 

1 

.9005 

. 1035xl0-5 

. 1035x10-5 

Mobile 

2 

.056 

a*  «• 

1 

.9477 

. 2464x10"2 

. 2464x10*2 

Las  Vegas 

808 

22.44' 

7.805 

38 

.9767 

. 1744x10-2 

.0663 

Belize 

2 

.056 

-- 

1 

.9477 

.8896x10-2 

.8896xl0“2 

Miles  City 

30 

.806 

.790 

2 

.9345 

1.5540x10-2 

3. 108x10-2 

San  Juan 

0 

-- 

-- 

0 

-- 

-  - 

0 

Column  1 
Column  2 
Column  3 

Column  4 

Column  5 
Column  6 
Column  7 


=  Total  number  of  gulps  in  36  years. 

=s  Mean  number  of  gulps  per  year. 

=  Standard  deviation  of  column  2.  Computed  only 
when  necessary 

=  Design  number  of  gulps  per  year.  Probability  ap¬ 
proximately  .  95  that  this  number  or  fewer  gulps 
will  occur. 

=  Computed  probability  of  design  number. 

=  Mean  water  gain  per  gulp,  grams  per  cubic  foot. 

=  Design  yearly  water  gain,  grams  per  cubic  foot. 
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CONCLUSIONS 


1.  The  modified  Monte  Carlo  method  used  in  this  study  provides  a 
quantitative  estimate  of  water  gains  by  containers  stored  in  the  open. 

2.  For  free  breathing  containers  tropical  and  monsoon  climates 
produce  the  greatest  water  gain,  with  a  mean  on  the  order  of  21  grams 
of  water  per  cubic  foot  per  year. 

3.  Of  the  U.  S.  stations  investigated,  the  largest  mean  gain  in  a  free 
breathing  container  was  14.  4  grams  per  cubic  foot  per  year  at  Mobile, 
Ala.  but  the  largest  single  year  gain,  31.2  grams  per  cubic  foot,  oc¬ 
curred  in  Las  Vegas,  Nev.  although  this  was  6  standard  deviations 
away  from  its  mean. 

4.  For  valved  containers,  greatest  gains  occur  in  a  desert  climate 
(Las  Vegas)  with  a  high  steppe  (Miles  City)  a  poor  second. 

5.  A  breather  charge  of  one  ounce  per  cubic  foot  should  provide  about 
six  months  life  in  world  wide  distribution  provided  the  container  and 
contents  are  adequately  dried. 

6.  All  valve  configurations  investigated,  including  the  +1/2,  -1/2  psig 
valve,  do  not  breathe  enough  water  to  affect  storage  life  of  the  con¬ 
tainer,  assuming  initial  desiccant  charge  conforms  to  MIL-P-116  and 
reasonable  storage  life  goals. 

7.  The  least  valving  configuration  which  would  not  adversely  affect 
storage  life  has  not  been  disclosed  by  this  study. 
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RECOMMENDATIONS 


1.  Free  breather  charge  be  set  at  one  ounce  per  cubic  foot  pending 
further  work  on  effect  of  breather  tubes  and  a  more  complete  determi¬ 
nation  of  the  distribution  function. 

2.  Work  be  done  on  breather  tubes  to  determine  their  effects.  This 
work  should  be  both  experimental  and  theoretical. 

3.  The  probability  distribution  function  of  water  gains,  for  three  or 
four  selected  stations,  be  determined  with  greater  precision  than  has 
been  possible  here. 

4.  Controlled  breathing  containers  be  accepted  with  valving  arrange¬ 
ments  as  low  as  +1/2,  -1/2  psig. 

5.  A  study  be  made  to  determine  what  the  least  pressure-vacuum 
spread  need  be. 

6.  Container  strength  requirements  be  based  on  the  pressures  and 
vacuums  developed  at  maximum  flow  rate  encountered  in  flight. 

7.  Maximum  rates  and  pressures  for  commercially  available  valves 
be  determined  and  incorporated  in  the  appropriate  container  design 


criteria. 
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APPENDIX  I 


DESICCANT  REQUIREMENTS 
FOR 

AIR  SHIPPED  CONTAINERS 


NOTE:  The  information  contained  in  this  Appendix  was  prepared  by 
M.  J.  Peterman  and  R.  S.  Nelson  of  Rocketdyne  Division  of  North 
American  Aviation.  It  was  furnished  to  Reed  Research  by  letter 
MJPfoc  dated  31  August,  1961. 

Gracious  permission  to  include  these  results  as  a  part  of 
the  overall  report  is  hereby  acknowledged. 


A.  INTRODUCTION 


The  AOA  container  design  section  Ad  Hoc  Committee  is  investigating 
pressure  and  desiccant  considerations  of  containers  for  air  trans¬ 
port  and  outside  storage.  Rocketdyne  was  assigned  the  subject  pro¬ 
ject.  Consequently  this  report  covers  theoretical  calculations  which 
answer  the  following  conditions: 

1.  How  much  water  per  cubic  foot  of  container  volume 
must  be  absorbed  by  the  desiccant  in  one  ascent  and 
descent  presuming  no  valve,  a  +2  -1  psi  valve,  or  a 
+1-1  valve? 

2.  What  is  the  relationship  between  valve  cracking 

pressures  and  water  influx? 

B.  SUMMARY 

Desiccant  requirements  were  calculated  for  containers  having  relief 
valves  with  various  settings  under  conditions  of  air  shipment.  Slow 
and  fast  rates  of  descent  were  considered. 

It  was  found  that  a  container  with  a  valve  having  a  three  psi 
difference  in  opening  and  closing  pressure,  and  under  conditions 
of  a  slow  descent,  would  require  approximately  ten  percent  less 
desiccant  than  an  unrestricted  breathing  container.  By  subject¬ 
ing  this  container  to  a  fast  descent  (2200  ft.  per  minute  used  as 
standard)  its  desiccant  requirements  would  increase  33%j  18% 
more  than  the  unrestricted  container  at  fast  descent.  To  verify 
these  findings,  calculations  were  made  of  a  container  having  a 
pin  hole  opening  (an  extreme  case  of  a  container  having  a  small 
opening  and  a  fast  rate  of  descent).  The  pin  hole  container  re¬ 
quired  49%  more  desiccant  than  the  unrestricted  free  breathing 
container. 

C,  CONCLUSIONS 

When  subjected  to  air  shipment,  a  container  having  a  relief  valve 
with  a  three  psi  differential  in  opening  and  closing  pressures, 
provides  little  more  moisture  protection  than  a  free  breathing 
container.  Small  changes  in  the  settings  of  the  relief  valves 
considered,  produce  no  significant  reduction  in  moisture  trans¬ 
mission  into  the  container. 
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The  closer  the  optimum  flow  pressure  of  the  valve  is  to  its  cracking 
pressure,  the  less  moisture  intake  into  the  container. 

The  difference  in  moisture  inhalation  between  a  container  with 
a  three  psi  differential  (+2,  -1)  and  one  with  a  two  psi  differen¬ 
tial  (+1,  -1)  is  quite  small  (5%). 

D.  DISCUSSION 


The  following  assumptions  were  used  as  the  basis  for  the 
calculations. 


1.  That  air  behaves  as  a  perfect  gas  according  to  Boyle's 
gas  law. 


2. 
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That  moisture  conditions  of  figure  IV,  of  MIL-STD-210A,  is 
a  valid  representation  of  the  actual  atmosphere  encountered. 


The  value  of  0.  022  pounds  of  moisture  per  pound  of  dry  air, 
representing  an  approximate  dewpoint  temperature  of  80°F 
is  taken  from  sea  level  to  the  altitude  (approximately  8,  000 
feet)  where  it  becomes  the  moisture  content  of  saturated  air. 
At  higher  altitudes  the  design  moisture  content  is  the  satura¬ 
tion  value  corresponding  to  the  temperature  and  pressure  for 
a  given  altitude.  Above  8,  000  feet  the  humidity  decreases 
quite  rapidly. 

3.  That  pressure  and  temperature  conditions  of  Chart  IV  of 
MLL-STD-210A  is  a  valid  representation  of  the  actual 
atmosphere  encountered.  This  chart  portrays  the  standard 
polar  atmosphere. 


4.  That  the  air  within  the  container  at  the  start  of  the  cal¬ 
culations  has  a  relative  humidity  of  40%. 

5.  That  during  slow  descent  conditions,  the  gage  pressure  of 
the  container  is  that  of  the  valve  cracking  pressure. 

6.  That  during  fast  descent  conditions,  the  gage  pressure  of 
the  container  is  that  of  the  valve  optimum  flow  pressure 
from  30,000  feet  to  sea  level.  At  sea  level  the  container 
pressure  rises  to  that  of  the  cracking  pressure. 
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7.  That  the  time  lag  in  water  adsorption  caused  by  desiccant 
adsorption  rates  is  not  a  significant  factor. 

it  is  realized  that  the  dessicant  will  not  absorb  the  water 
instantaneously  and  some  lag  will  occur.  It  is  felt  that  this 
lag  will  be  relatively  short  and  of  little  importance. 

Due  to  the  *a~t  that  data  necessary  for  the  calculations  was 
in  the  form  of  graphs,  it  was  impossible  to  calculate  directly 
the  influx  of  water  by  means  of  mathematical  equations. 

The  amount  of  water  influx  wao  determined  by  averaging  the 
pvassnze,  temperature,  and  humidity  changes  in  3,000  foot 
increments,  Largev  increments  could  not  be  used  because  of 
the  exponential  shape  of  the  curves. 

''he  method  of  ca] emulation  waa  as  follows: 

1.  The  atmospheric  pressure  on  the  outside  of  the  container  was 
decreased  to  that  of  30,  000  feet.  The  inside  pressure  was 
calculated  from  the  valve  cracking  pressure.  Assuming  that 
the  container  temperature  drops  to  that  of  the  outside  air,* 
the  pounds  of  air  in  the  container  were  then  calculated. 

*Air  temperature  within  the  plane  is  greater,  although 
unknown,  resulting  in  an  additional  safety  factor. 

2.  The  atmospheric  pressure  on.  the  outside  of  the  container  was 
raised  in  3,  000  foot  increments  until  the  valve  opened  and 
admitted  car.  The  resultant  temperature  inside  the  container 
was  calculated  by  assuming  adiabatic  mixing  with  the  incoming 
air  and  in  proportion  to  the  pressure  change  inside  the  con¬ 
tainer. 

3  The  pounds  of  air  at  the  end  of  each  increment  were  cf  iculated 
using  the  resultant  temperature  and  pressure.  From  thic,  tne 
pounds  of  air  entering  the  container  were  derived, 

4,  The  average  altitude  fer  eucn  increment  was  calculated,  and 
with  this,  the  humidity  for  the  increment  obtained  from  Figure 
IV  of  MIL-STD-210A. 

5.  The  pounds  of  water  picked  up  in  each  increment  was  calculated 
by  multiplying  the  aveiage  humidity  by  the  pounds  of  air  in¬ 
haled. 


1-4 


6.  The  total  moisture  to  be  adsorbed  by  the  desiccant  was  ob- 
tained  by  adding  the  moisture  inhaled  in  each  increment. 

7.  The  pounds  of  moisture  was  then  converted  to  grams.  The 
units  of  desiccant  required  was  calculated  by  dividing  the 
grams  of  water  by  six.  To  maintain  a  relative  humidity  of 
40%  each  unit  of  desiccant  (MIL~D»3464)  must  absorb  6 
grams  of  water. 

E.  SUMMARY  OF  RESULTS 

A  1,000  cubic  foot  container  requires  the  following  amounts  of 
desiccant  for  one  cycle  of  ascent  and  descent  from  30, 000  feet 
sea  level  and  maintain  a  relative  humidity  of  40%  within  the  container. 


Free  Breathing  Container 
(Very  large  opening) 

Free  Breathing  Container 
(Pin  hole  opening) 

Container  with  Valve 
(Cracking  pressure  +1,  -1  PSI) 

Container  with  Valve 
(Cracking  pressure  42,  -1  PSI) 

♦(Optimum  Flow  Pressure  +4.  5,  -2.  5  PSI) 

Container  with  Valve 
(Cracking  pressure  +3,  -0  PSI) 


Slow  Descent  Fast  Descent 

(2200  feet  per  minute) 

71  units  71  units 

71  units  106  units 

67  units  -*■» 

64  units  84  units 

62  units 


♦The  optimum  flow  pressure  is  the  required  pressure  differential  across 
the  valve  necessary  to  allow  its  rated  flow  volume. 
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